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Abstract: We study the impact of parton correlations on the double Drell-Yan process, 
i.e. on the production of two electroweak gauge bosons by double parton scattering in a 
single proton-proton collision. Spin correlations between two partons in a proton are shown 
to change the overall rate of the process and to induce characteristic angular correlations 
between the decay leptons of the two gauge bosons. 
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1 Introduction 

In the era of the LHC, a thorough understanding of the strong-interaction dynamics in 
proton-proton collisions is an imperative. Whereas the determination of parton distribution 
functions and the computation of hard-scattering processes at parton level are becoming 
fields of precision physics, other aspects remain much less well understood, both at a 
conceptual and a practical level. Among these aspects are multiparton interactions, where 
in a single proton-proton collision more than one parton in each proton participates in a 
hard scattering. Multiparton interactions can give sizeable contributions both to signal 
and to background processes, notably to Higgs production [1], to electroweak processes 
[2-10] and to multijet production [11-18]. They are also of interest in their own right from 
the point of view of hadron structure, because they contain information about correlations 
between partons inside the proton. 

Experimental evidence for hard multiparton interactions has been found at the ISR 
[19], the SPS [20] and the Tevatron [21-25]. Due to the rapid increase of parton densities at 
small momentum fractions, one expects these interactions to be even more prominent at the 
LHC, and first results have indeed been reported [26], with more studies being expected 
in the future [27]. Substantial effort has gone into the modeling and implementation 
of multiparton interactions in Monte-Carlo generators [28-34]. A mini-review on earlier 
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developments of the subject is given in [31], and an overview of current developments can 
be found in the conference proceedings [35, 36]. Recently there has been renewed interest in 
understanding the theoretical foundations of multiple hard scattering [37-47] , with many 
issues remaining to be clarified or worked out. 

The simplest assumption one can make when computing multiple hard scattering is 
that the different partons in the proton are uncorrelated with each other. This leads to 
very compact results, but it certainly is a simplification whose validity and limitations 
need to be investigated. Recent phenomenological studies of different types of correlation 
effects in multiparton interactions can be found in [17, 34, 48-52]. In this work we focus on 
correlations between the polarization of two partons in an unpolarized proton and on their 
consequences for the overall rate and for final-state distributions in double hard scattering 
processes. We build on the observations made in [40] and extend the results of [46]. The 
relevance of spin correlations for multiple interactions was pointed out long ago in [53] , but 
that work have not been followed up until recently. 

We perform our study for double Drell-Yan production, where two electroweak gauge 
bosons (j* , Z,W^) are produced in two independent quark-antiquark annihilation pro- 
cesses. This has long been recognized as a prototype for multi-parton interactions [54-56]. 
Unfortunately, recent phenomenological estimates [8, 9, 57] find that rates for double Drell- 
Yan production at the LHC will likely be too small to allow for a detailed experimental 
investigation of final-state distributions. We nevertheless choose this process for our inves- 
tigation, given that it is among the simplest double scattering process from the theoretical 
point of view but still exhibits a wealth of nontrivial features. The results we obtain for 
double Drell-Yan production will have their analogs for processes with higher rates, such 
as the production of four jets or the production of one gauge boson plus a dijet. (Note 
that up to a global factor the tree-level graph for the production and subsequent decay of 
a gauge boson is identical to the graph for qiq2 —> </3<?4 with a gluon in the s-channel.) 

Even if two partons in a proton are unpolarized, there can be correlations between 
their transverse positions and their longitudinal momentum fractions. There are indeed 
reasons to expect such correlations, see e.g. [34] and section 2.6 of [40], but not much 
is currently known about them. In the spirit of an exploratory study, we will calculate 
these correlations in a simple model of the proton and discuss their consequence for the 
transverse-momentum spectrum of double Drell-Yan production. 

This paper is organized as follows. In the next section we recall the basic formalism 
for computing double parton interactions, define double parton distributions that describe 
spin correlations between the two partons and give the polarization dependent parton-level 
cross sections. In section 3 we give our results for the cross section of the double Drell-Yan 
process. In section 4 we investigate the correlation between transverse and longitudinal 
variables in double parton distributions within a simple model. Section 5 summarizes our 
findings, and in an appendix we list the coupling factors entering the cross section formulae 
in section 3. 



- 2 - 



2 Double parton scattering 



Consider the production of two gauge bosons Vx, V2 = 7*, Z, W in a pp collision, followed 
by the leptonic decays j*,Z — > £ + £~ or W — > £v. Four- momenta are assigned as p(p) + 
Pip) ~^ Viill) + ^2(92) + X. We are interested in the fully differential cross section of the 
four-lepton final state and restrict ourselves to the kinematic region where the transverse 
momenta q-^ and q 2 of the gauge bosons in the pp center-of-mass are much smaller than 
their invariant masses, i.e. we assume q 2 ,q 2 <C <?i,#2- ^ ^ s i n this kinematics that double 
parton scattering is not power suppressed compared with the production of the gauge 
boson pair by a single hard scattering [39]. In the calculation of the cross section, the 
invariant masses Qi = iqf) 1 ^ 2 will serve as the hard scale necessary for the application of 
factorization. For simplicity we shall not assume any particular hierarchy in size between 
q 1 and q 2 or between Q\ and Q 2 . 

We assume that the double hard scattering cross section factorizes into the product 
of a double parton distribution (DPD) in each proton and a parton-level cross section for 
each of the two hard scatters. This factorization has not been proven, but several elements 
of such a proof have been given in [40] . Schematically, the double parton scattering cross 
section then reads [39, 40] 



da 



Yli-i dxi dxi d 2 qi dQi 



= 1 da! da 2 f d l z x d 2 z 2 iziqi - iZ2q2 
DPS C dUi dU 2 J (2tt) 2 (2tt) 2 

x J d 2 y F(x 1 ,x 2 ,z 1 ,z 2 ,y) F(x 1 ,x 2 ,z 1 ,z 2 ,y) (2.1) 



with a combinatorial factor C equal to 2 when the final states of the two hard interactions 
are identical and equal to 1 otherwise. Here ddijd^li denotes the cross section for quark- 
antiquark annihilation into a lepton pair via the gauge boson Vi, taken differential w.r.t. 
the lepton angles in the appropriate boson rest frame (see section 2.2). In the pp center- 
of-mass we define the z axis to point into the direction of the proton momentum p and 
use light-cone coordinates = (v° db u 3 )/\/2 and v = (u 1 ,^ 2 ) for any four- vector v. 
The kinematic variables x i = qf /p + and x i = q^/p~ determine the longitudinal parton 
momentum fractions in the DPDs, which we denote by F for the proton with momentum 
p and by F for the proton with momentum p. The arguments Zj and y of the distributions 
determine where the hard-scattering processes take place in transverse configuration space. 
As indicated in figure 1, y is the transverse distance between the two scattering partons in a 
proton (and hence between the two annihilation processes) if one takes the average between 
the scattering amplitude and its conjugate. As shown in [39, 40] z% is the Fourier conjugate 
variable to the transverse momentum of parton % (again averaged between amplitude and 
conjugate amplitude). F(xi, x 2 , z\, z 2 ,y) is thus the Fourier transform of a transverse- 
momentum dependent DPD. The factorization formula (2.1) generalizes the expression 
for single Drell-Yan production in terms of transverse-momentum dependent single-parton 
densities [58, 59]. 



Figure 1. A graph for the double Drell-Yan process, where two quarks in the right-moving pro- 
ton interact with two antiquarks in the left-moving proton. The figure shows the assignment of 
four-momenta (p, p, qi, 52), of light-cone momentum fractions (xi, Xi) and of transverse position 
arguments as explained in the text. The dotted vertical line denotes the final-state cut. 

Integrating (2.1) over the transverse boson momenta q-y and q 2 one obtains 

Here F(xi,X2,y) and F(x\,x<2,y) are transverse-momentum integrated (also called col- 
linear) DPDs, which were introduced long ago in [55, 60]. Naively, they are obtained by 
setting z\ = 22 = in the distributions that appear in (2.1). Closer analysis reveals that 
transverse-momentum dependent and collinear DPDs (just as their counterparts for single 
partons) require different types of regularization and subtractions of divergences. As a 
result the distributions depend in different ways on an ultraviolet renormalization scale, 
and (with the exception of specific distributions) also on a rapidity parameter which is 
closely related to Sudakov logarithms. This is discussed in [40] and will be tacitly implied 
in the remainder of the present work. 

Equations (2.1) and (2.2) are schematic in that they omit labels for and summation over 
the quantum numbers of the partons (quarks vs. antiquarks, flavor, polarization and color). 
This information will be restored in section 2.3. We emphasize that these equations only 
give one contribution to the cross section for four-lepton production. Further contributions 
need to be added from the familiar single hard-scattering mechanism (where the four 
leptons are produced in a single parton-level process), the interference between single and 
double hard scattering, as well as double hard-scattering graphs with fermion number 
interference [40]. The single hard-scattering contribution is straightforward to compute 
(see e.g. [57]), whereas the different interference contributions are not. As argued in [40], 
the fermion number interference contribution should become relatively unimportant at 
small momentum fractions X{, Xj. 



da 



IIL1 dxi dxi dVLi 
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2.1 Double parton distributions 

We now briefly specify the DPDs necessary for our study, referring to [40] for more detail. 
For two quarks in an unpolarized right-moving proton we write 

F aia2 (x u x 2 , Zl ,z 2 , y) = 2p + j ^ dy~ e ^P + +^r + {p \ o a , 2 (0, z 2 ) O ai (y, z x ) \p) , 

(2.3) 

where an average over the proton polarization is implied. The operators 

Oa i {y,z i ) = q i {y-\z i )T ai q i {y^\z i ) (2.4) 

z, T=y =o 

for quarks of flavor qi are understood to include appropriate Wilson lines between the two 
quark fields, as well as appropriate regularization and subtractions as mentioned above. 
Notice the correspondence of position arguments in (2.3) and (2.4) with figure 1. The Dirac 
matrices 

r, = ±7 + > r Ag = i 7 +7 5 , r^ = i^ + 75 (j = 1,2) (2.5) 

select different polarizations of the quarks in the proton, with q corresponding to unpolar- 
ized quarks, Aq corresponding to longitudinal polarization and 5q to polarization in the 
transverse direction j. Notice that the labels in (2.3) specify both the flavor and the 
polarization of the quarks. In full analogy one can define DPDs F^ o 2 f° r two antiquarks, 
as well as quark-antiquark distributions i*ki,a 2 an d F ai ^ 2 - 

The quark coupling to gauge boson V\ need not have the same flavor in the scattering 
amplitude and its complex conjugate, because a mismatch in flavor can be compensated by 
the quark coupling to gauge boson V 2 . For example, the quarks with transverse positions 
y + \z\ and —\z 2 in figure 1 can be u quarks if the quarks with transverse positions \z 2 
and y — \z\ are d quarks. The DPDs describing this type of quark flavor interference are 
given by 

Fl ia2 { Xll x 2 ,z llZ2l y) = 2p + J d -h- d -^-dy~ e *w + +- 2 ^>+ ( p | Oi 2 (0, z 2 ) J ai ( y , Zl ) \p) 

(2.6) 

with the product of operators 

O I a2 (0,z 2 )O I (y,z 1 )=q 1 (-lz 2 )T a2 q 2 (lz 2 )q 2 ( y -^z 1 )T ai q 1 ( y + ^z 1 ) . 

z i = 2 2 =y =° 

(2.7) 

We note that these distributions are complex valued and that their imaginary part changes 
sign when one interchanges the flavor (but not the spin) assignments and replaces Zj — > —Zi, 
e.g. 

F q 1 Aq 2 ( x l,X2,z 1 ,z 2 ,y) = [F^ Aqi (x 1 ,x 2 ,-z 1 ,-Z2,y))* ■ (2.8) 
As we shall see, this ensures that physical cross sections are real-valued. 
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Let us now classify the different combinations of quark polarization, taking into account 
the constraints of parity invariance [40]. For unpolarized and longitudinally polarized 
quarks we have 



F qg = f qq (xi,x 2 ,zi,z 2 ,y) , 

FqAq = g q Aq(xi,X2,Zi,Z 2 ,y) , 



FA q A q = fAqAq{xi,X 2 ,Z 1 ,Z 2 ,y) 
FAqq = 9Aq q (xi,X 2 ,Zl,Z 2 ,y) , 



(2.9) 



where / denotes scalar and g pseudoscalar functions. For transverse quark polarization the 
parton distributions carry a transverse index and can be decomposed as 



where the scalar and pseudoscalar functions depend on the same variables as in (2.9). Here 
y % = e %3 y 3 is a transverse vector orthogonal to y l , defined in terms of the two-dimensional 
antisymmetric tensor e %3 (with e 12 = 1). Factors of the proton mass M have been introduced 
in order to have the same mass dimension for all distributions / and g. For two transversely 
polarized quarks we finally write 



Decompositions analogous to (2.9) to (2.11) hold for antiquarks and for flavor interference 
distributions. 

Corresponding definitions apply for two partons in a left-moving proton, with + and — 
components interchanged in (2.3) to (2.7). Note that the covariant expression of the two- 
dimensional antisymmetric tensor in terms of the four-dimensional one is e 13 = e" 1 13 (with 
e oi23 = !)• bi the analogs of (2.9) to (2.11) for left-moving partons one hence needs to 
change the sign of y and of the pseudoscalar functions g (which can be written as e 13 
contracted with a parity even tensor constructed from z±, z 2 and y). 

All distributions discussed so far allow for two color structures, one where the two 
fields in the operator O ai are coupled to a color singlet and one where they are coupled to 
a color octet. This requires a further index on all distributions, which we will not display 
in the present work for brevity. 

2.2 Reference frames 

Let us now introduce the reference frames and coordinate axes needed to describe the 
angular dependence of the cross section. 

In the pp center-of-mass we have the z axis pointing along the momentum p. The 
four- vector defining this axis is hence = (p — p)^/\/2pp, where here and in the following 
we neglect the proton mass. We choose a fixed four- vector orthogonal to p and p to 
define the x axis. The precise choice does not matter for our purpose, but one may for 
instance adopt the convention to have the x direction point towards the center of the LHC 





(2.10) 



F7 qSq = S ij f SqS q + M 2 (2yV " y 2 ^) f 5 qSq 

+ M 2 {tfy 3 + tfy 3 ) g s 5q5q + M 2 {y>y 3 - tfy 3 ) g a SqSq . 
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k 



lepton plane 



Figure 2. Coordinate system in the rest frame of vector boson Vi. The z axis bisects the angle 
between the spatial components of the momenta p and — p, and the x axis corresponds to a fixed 
reference direction as explained in the text. (In general, the proton momenta are therefore not in 
the x-z plane.) U and U are the momenta of the lepton and the antilcpton from the boson decay, 
respectively. 9i denotes the polar and ifi the azimuthal angle of the lepton. Note that Lpi is negative 
in this example. 

ring. The y axis is then denned such as to obtain a right-handed coordinate system; the 
corresponding four-vector can t>6 written as — vpo 

X v pPp a I (pp). 

The kinematics of the gauge boson decays into lepton pairs is conveniently described 
in the rest frame of the respective gauge boson. The z axis in the rest frame of the boson 
Vi is defined by the four- vector 

% = \\/<% + *i 

where is the transverse boson momentum in the pp center-of-mass as before. As illus- 
trated for one boson in figure 2, the z axis bisects the angle between the spatial components 
of p and — p in the boson rest frame. The x axis is specified by 



and the y axis is again defined so as to obtain a right-handed coordinate system, i.e. by 
= e^vpo- X? qf/Qi- With these reference axes we define the polar and azimuthal 

angles Qi and ifi of the lepton (as opposed to the antilepton) in the decay of Vi, i.e. of t~ 

in the decay of a 7* , Z or W~ and of V(_ in the decay of a W + . 

Noting that Xqi is the x component of qi in the pp center-of-mass, we see in (2.13) 

that X\, X2 and X differ from each other by amounts of order {q^/Q^ which is a small 

parameter in our calculation of the cross section. Likewise, one finds differences of order 

\li\/Qi between Y±, Y% and Y. As we shall see shortly, this greatly simplifies the discussion 

of azimuthal angles in our calculation. 

Readers familiar with the analysis of single Drell-Yan production will recognize that 

our choice of z axes is the same as in the Collins-Soper frame [61]. Useful information 
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p_ 



(2.12) 



X» 



Xqj 

Qi 



(2.13) 
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about this frame can e.g. be found in [62, 63]. By contrast, we define x axes (and thus the 
azimuthal angles ipf) starting from a fixed direction in space, whereas in the Collins-Soper 
frame the x axis is defined such that the proton momenta lie in the x-z plane. The latter 
choice becomes undefined when the transverse boson momentum in the pp center-of-mass 
goes to zero. For unpolarized single Drell-Yan production this is not a problem because in 
this limiting case all azimuthal dependence in the cross section must vanish due to rotation 
invariance. However, in the double Drell-Yan process there can be an azimuthal dependence 
even if q l or q 2 or both go to zero, as we shall see. Choosing one of these vectors (or any 
linear combination of them) to define x axes would therefore entail ill-defined azimuthal 
angles at some point in phase space where there can be a nontrivial azimuthal dependence. 

The physical cross section must of course not depend on the arbitrary fixed direction 
specified by X^ 1 . To understand how this happens, we anticipate that our results will 
depend only on the difference of azimuthal angles whose definition depends on X^, such 
as for instance ipi — ipi- These angles are defined in different frames, but to the accuracy 
of our calculation we can replace them with the azimuthal angles of the leptons in the pp 
center-of-mass. This is easily seen by writing trigonometric functions of ipi in terms of 
invariant products X^li and Y^Zj, where li is the four-momentum of the lepton from the 
decay of V{. When calculating the cross section we neglect terms of order \qi\/Q i and can 
hence approximate Xili ~ Xli and Y%1% ~ Yli, which gives azimuthal angles in the pp 
center-of-mass as announced. 

2.3 Hard-scattering cross sections 

If we restore labels for quarks and antiquaries, their flavor and their polarization, the cross 
section formula (2.1) reads 

da 1 v-^ f d 2 z\ d 2 z 



Y\ i=1 dxidxid 2 qidQi C ^ 



C L J (2n) 2 (2n) 2 ] dV 



a 2 a 3 a 4 



da aia3 da a2a4 - ^ da aia3 da a2 g 4 - 

£ aic^-f a3<i4 T ,„ ,„ " aia2 r aACli 



dni dn 2 1 dtti dti 2 

da aias da a2 a4 jp f? * da^as daa 2 a4 



+ {flavor interference} , (2-14) 

where here and in the following we omit the label "DPS" for double parton scattering. 
In all terms the DPDs have arguments as in (2.1), which will be omitted henceforth for 
brevity. To distinguish the distributions for the left- and right-moving proton we use the 
notation F and F and a corresponding notation for the scalar and pseudoscalar functions 
/, / and g,g introduced in section 2.1. The first subscript in da aiaj and daa iaj denotes the 
right-moving parton and the second subscript the left-moving one. The sum over a\ to 04 
runs over quark flavors and polarizations (q, Aq,5q). 

The flavor interference terms involve the interference DPDs in (2.6) and corresponding 
interference terms for the hard scattering. These interference terms only appear if the 
produced bosons are both neutral or both charged, otherwise the quark and antiquark 
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flavors in the annihilation processes do not match. We will return to this in the next 
section. 

Labels for the color structure of the DPDs are not displayed in (2.14). With the con- 
ventions of [40], each factor F F is to be replaced with the sum 1 F 1 F + 8 F 8 F of color singlet 
and color octet distributions, without change in the hard-scattering cross sections. This 
holds for the production of arbitrary color-neutral states in the hard-scattering processes. 

It is straightforward to compute the tree-level cross section for quark-antiquark anni- 
hilation into a gauge boson followed by its leptonic decay. In accordance with the power 
counting scheme underlying the cross section formula (2.14), the transverse boson momenta 
q i are set to zero in this calculation since by assumption they are small compared with 
the invariant mass Qi. This also simplifies the kinematics of the gauge boson decays as we 
already noticed in section 2.2. 

Consider first the case where both quark and antiquark are unpolarized or longitudi- 
nally polarized. The angular dependence of the cross section is then of the form 

= (l + cos 2 ^)^.^) +2coB0 i K' a .a s (Qi), (2.15) 

with di = qi,Aqi and dj = q~j,Aq~j. The integration element reads dfli = difidcosOi as 
usual. The factors K and K' depend on coupling constants and on Qi via the gauge boson 
propagators. One easily finds 

KAqiAqj = ~K qi q. , K q % Aq, = ~^Aq iqj (2.16) 

and analogous relations for K', so that 

da qi q. _ da Aqi A qj da qi A qj _ daA qiqj 



dVli dQj d£lj dttj 



(2.17) 



Because of chirality conservation for massless quarks one has vanishing parton-level cross 
sections for the annihilation of a transversely polarized parton with an unpolarized or 
longitudinally polarized one, das qiqj = das qi Aq 3 = da qi s qj = daA qi Sqj = 0. If both quark 
and antiquark are transversely polarized, one finds 



d(T S l q t Sqj ■ 2 

— sin 



dUi 

+ 



|[cos(2^) K SqiS q 3 (Qi) - sin(2^) K' 5qi5q . (Q t ) ] (X k X l - Y k Y l ) 

sin(2^) K 5qiSqj (Qi) + cos(2 W ) K's^iQi)] {X k Y l + Y k X 1 )} (2.18) 



with X and Y as defined in section 2.2. The transverse indices k,l in (2.18) refer to the 
pp center-of-mass, where they are to be contracted with the corresponding indices of the 
DPDs. We note that contraction of (2.18) with the transverse spin vectors s k , s l of the 
quark and the antiquark gives the simple expression 

da kl 

d Q qj s k s l = sin 2 Qi [cos(<^ s + tp s - 2tpi) K 5qi s q . + sin(v? s + ip s - 2<pi) K' Sq . Sq .] , (2.19) 
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where (p s and <p§ are the azimuthal angles of the spin vectors in the pp center-of-mass and 
our normalization convention is s 2 = s 2 = 1. 

The preceding expressions hold for both neutral and charged vector bosons, and the 
coupling factors and K' ai3ij appearing in (2.15) and (2.18) are given in appendix A. 

For neutral boson production the annihilating quark and antiquark have the same flavor. 
In this case we will use da qiq . , da qiAqj , . . . with i 7^ j to denote the interference terms for 
flavor qi in the amplitude and flavor qj in the conjugate amplitude. The relations (2.15) 
to (2.19) remain valid for these interference terms. As can be seen in appendix A, the 
corresponding coupling factors are complex, and their imaginary parts change sign when 
the flavor (but not the spin) labels are interchanged, e.g. 

^<J192 = (-^92 91 )* j -^9iA92 = (Kq 2 Aqi)* • (2.20) 

We note that for invariant masses Qi far below the Z mass, the neutral boson channel is 
well approximated by 7* production alone. The only nonzero coupling factors in this case 
are K qiqj = -K AqiAq . = K SqiSqj . 

For W boson production we use da qiqj , da mAqj etc. to denote cross sections with 
different flavors qi, qj in the initial state. We do not need a separate notation for flavor 
interference terms in this case, because the product da ai a 3 dd a2 a A of cross sections for WW 
production is equal to the product of the corresponding interference terms, except for 
CKM factors that can easily be identified. Using that W bosons only couple to left-handed 
fermions, we find further simplifications for the coupling factors: 

Kqtqj = KqiAqj , ^SqiSqj = 0, ^oia 2 = K ai a 2 , (2-21) 

where the second relation reflects that the operator 0$ q for transverse quark polarization 
corresponds to the interference between left- and right-handed quarks. Together with the 
relations (2.16) we are thus left with only one independent coupling factor for W~ and 
only one for W + production. 

So far we have discussed cross sections and interference terms da ai aj for the annihilation 
of a right-moving quark with a left-moving antiquark. The cross sections and interference 
terms daa jai f° r right-moving antiquarks and left-moving quarks have the same form as in 
(2.15) and (2.18). The associated coupling factors are given by 

K- mi = (K qiqj y , K' qjqi = -{K' qiqj )* (2.22) 

and analogous relations for the spin combinations At/Ag and Sq5q, and by 

KqjAqi = - {K qi Aqj)* , K 'q 3 A qi = {K' qi Aq 3 ) * (2-23) 

and an analogous relation for the spin combination Aqq. 

3 The double Drell-Yan cross section 

Inserting the hard-scattering cross sections (2.15), (2.18) and the DPD decompositions 
(2.9) to (2.11) into the factorization formula (2.14), we obtain our final results for the 
double parton scattering contribution to four-lepton production. 
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For the production and decay of two W bosons, the result has a simple structure thanks 
to the relations (2.21), 



da 



WW 



Ui^ l dx i dx i d 2 q i dn i C ^ 
x (1 + cos0i) 2 (l + cos6> 2 ) 2 

X \fqiq2 + fAq\Aq 2 ~~ 9q\Aq 2 ~ 9Aqiq 2 ) ifq 3 q 4 + fAq 3 Aq 4 ~ 9q 3 Aq 4 ~ 9Aq 3 q 4 ) 
+ (1 + COS(9l) 2 (1 - COS 9 2 f 

X \fqiq~4 ~ /AqiAq4 + 9q\Aq 4 ~ 9Aqiq 4 ){fq 3 q 2 ~ fAq 3 Aq 2 + 9q 3 Aq 2 ~ 9Aq 3 q 2 ) 

+ (1-cos6 1 ) 2 (1 + cos9 2 ) 2 

X (/<?392 ~ fAq 3 Aq 2 ~ 9q 3 Aq 2 + 9Aq 3 q 2 ){fqiq 4 ~ /AgiAg 4 ~~ 9q 1 Aq 4 + 9Aqiq 4 ) 

+ (l-cos6ii) 2 (l-cos6>2) 2 

x (/g 3 f 4 + fAq 3 Aq 4 + <7g 3 Ag 4 + 9Aq 3 q 4 ){f qi q 2 + fA qi Aq 2 + 5giA<j 2 + 9A qi q 2 ) 

+ {flavor interference} , (3-1) 

where the sum over g 4 to 54 runs over quark flavors. The flavor interference terms are 
obtained by replacing the DPDs in one or in both protons with their interference analogs 
and by appropriately changing the CKM factors in the product K qiq3 K q2qi . Different types 
of flavor interference terms are shown in figure 3. Taking into account the minus sign in the 
definition of pseudoscalar distributions for left-moving partons, e.g. in F gi Ag j = —g qi A qj , 
we recognize that the DPD combinations in (3.1) correspond to negative-helicity quarks 
and positive- helicity antiquarks, as required by the left-handed nature of the charged weak 
current. 

We see that for W pair production the presence of longitudinal parton spin correlations 
in the proton changes the overall rate of the cross section as well as the distribution in the 
polar angles of the decay leptons. 

For one or two neutral bosons (7*, Z) the structure of the cross section is more com- 
plicated. We split the cross section (2.14) into three parts, for the case without trans- 
verse quark polarization and for the cases where one or two hard interactions 
are initiated by transversely polarized quarks. The contribution with only unpolarized and 
longitudinally polarized partons reads 

1 ^ f d 2 Zl d 2 z 2 



C „hJ (2^) 2 (2vr) 2 J V 



Y\ i=1 dx i dx i d 2 q i dn i C q 

. Tf' _ 1 171 l ™c,2 Q„\TS _ I O^ofl- I. 

V <?2 94 



{ [(1 + cos 2 9 1 )K qi g. A + 2 cos #1 K' qi g 3 ] [(1 + COS 2 02)K q2 g 4 + 2 COS 6 2 K' q2 

X (/<?l 92/93 94 + fAqiAq 2 fAq 3 Aq 4 + 9qiAq 2 9q 3 Aq 4 + 9Aqiq 2 9Aq 3 q 4 ) 
+ [(1 + COS 2 0{)K qi Aq 3 + 2 COS 9 1 K' qiA q 3 ] [(1 + COS 2 2 )K q2 Aq 4 + 2 COS 9 2 K' q2 Aq 4 
X {fqiq2fAq 3 Aq 4 + f'Aq 1 Aq 2 fq 3 q 4 + 9q\Aq 2 9Aq 3 q 4 + 9Aq 1 q 2 9q 3 Aq 4 ) 
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Figure 3. Hard-scattering graphs for the production of W + W + (a, b) or of W + W~ (c. d). The 
labels q and q indicate whether a parton corresponds to a quark field or a conjugate quark field in 
the relevant DPD. Graphs (b) and (d) are multiplied with interference distributions for one of the 
protons, whereas graphs (a) and (c) go along with interference distributions for both protons. 

- [(1 + cos 2 9i)K qi g a + 2 cos 0i K' qiqa ] [(1 + cos 2 6 2 )K q2Aqi + 2 cos 9 2 X£,A&] 

x \9q1Aq2fq3q4 + QAq^fAq^Aqi + fq^QqzAqA + /Agi AiJ2<? A5354) 

- [(1 + cos 2 e^Kq^ + 2 cos 6 1 K' qiA q 3 ) [(1 + cos 2 d 2 )K q2 g 4 + 2 cos 6 2 K' q2q4 ] 

x {9qiAq 2 fAq 3 Aq4 + 9Aq 1 q 2 fq 3 q4 + fqiq29Aq 3 qi + /a<ji Aq 2 9q^ Ag 4 ) | 

+ {flavor interference} + {gg permutations} . (3-2) 

The qq permutation terms are obtained by permutation of the quark-antiquark assignments 
in the DPDs and in the coupling factors K, K' as specified in (2.14). For neutral bosons 
the annihilating quark and antiquark have the same flavor, i.e. one has q\ = q% (q 2 = 94) if 
Vi (V2) is neutral. The flavor interference term for neutral boson pairs is then obtained by 
replacing all distributions f,g,f,g with their interference analogs / , g , f , g and by 
interchanging 1 f> 2 in the second subscript of the coupling factors, e.g. K qiqi K' q2 /\ q2 — > 
K qiq2 K q2 /± qi . The relations (2.8), (2.20) and their analogs for other polarizations ensure 
that the sum over all flavor assignments in (3.2) gives a real-valued cross section. Example 
graphs for flavor interference are shown in figure 4. 

We see in (3.2) that longitudinal parton spin correlations change the overall rate of 
double parton scattering and the dependence on the polar angles of the leptons, due to 
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Figure 4. Hard-scattering graphs for the production of two neutral gauge bosons. The labels q 
and q have the same meaning as in figure 3. 

the differences between the coupling factors K qiq ., K' qiqj and K qi ^ qj , K' qi & q .. Only in the 
neutral boson channel at Qi values small enough to neglect Z production does one have a 
fixed angular dependence da^ jd cos 0j oc 1 + cos 2 0j. 

We now turn towards the part of the cross section where one of the two annihilation 
processes involves transverse quark polarization (and thus produces a neutral gauge boson). 
It reads 

UU dxi dxi d* qi dQ t C ^ J (2vr) 2 (2vr) 2 J 

x ( [(1 + cos 2 O^Kg^ + 2 cos 0i i^ l9 - 3 ] 

X | [ COS 2(992 - fy) K Sq2 Sq 2 ~ Sm 2{lf 2 ~ fy) K6q 2 5q 2 ] 

X ifqi8q2 fqz&q~2 ~ 9qi8q29qs8q2 ~ f Aqi5q2 f Aqs5q2 9 Aqi&q29 Aq-iSq 2 ) 
+ [sin 2(992 - fy) Ks q2 Sq 2 + COS 2^2 - fy) K's q2 Sq 2 ] 

X ifqiSq29qa8q2 9q\8q2fqz&q~2 J 'AqiSq29 'Aqz5q~2 9Aq\5q2 J Aqz5q~2 ) ^ 
- [(1 + COS 2 9l)K qiAq3 + 2 COS 01 KgjAgs] 

x j [ cos 2(992 - 99^) i^ q2 <5q 2 - sin 2(992 - f y ) K' Sq2 s q2 ] 

X {fqiSq29Aq~38q~2 ~ 9q\bq2 f Aq~3&q~2 ~ f Aqi8q29q~3,&q~2 9Aq\8q2 /g3<5<?2 ) 
+ [sin 2(992 - fy) K Sq2 8q 2 + COS 2(992 - fy) K's q2 S q2 ] 

X (fqi6q2fAq~38q2 + 9q\8q29 Aq^q2 + f Aq\&q2 Iq-iSq2 + 9 Aqi5q29q~35q2 ) ^ ) 

+ {flavor interference} + {qq permutations} + {transv. pol. in interaction 1} , (3.3) 

where the flavor interference and qq permutation terms are obtained in the same way as in 
(3.2). The terms for transverse polarization in interaction 1 are obtained by replacing labels 
asl— t-2, 2— t-1,3— t-4 in the coupling factors and by making the same replacement in the 
DPD subscripts after interchanging their order, i.e. f aiSq2 -)• fs qia2 , 9a 3 5q 2 ~> 5%a 4 etc. 
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The azimuthal angle ip 2 of the lepton produced in interaction 2 has already been de- 
fined, and ipy is the azimuthal angle of y in the pp center-of-mass. As anticipated in 
section 2.2, the cross section depends only on the difference ip 2 — fy of these angles, in 
agreement with rotation invariance. The tp y dependence in (3.3) arises from the uncon- 
tracted vectors y and y in the DPDs (2.10) for transversely polarized partons: it is hence 
this polarization which enables a dependence of the cross section on the azimuthal angle 
of the produced lepton. 

The transverse distance y is integrated over in (3.3) and hence not measurable. The 
y integration is nontrivial because the DPDs depend on the azimuthal angles between y 
and Z\ and z 2 , whose directions are in turn correlated with those of q^ and q 2 through the 
exponential eT %Ziqi ~ lZ2tl2 . ^he integral over y, z\ and z 2 m the cross section thus turns the 
(fy dependence into a dependence on the azimuthal angles of the transverse momenta q^ and 
q 2 . All together we thus see that a correlation between y and the transverse polarization of 
parton 2 in the DPDs leads to an azimuthal correlation between the lepton from interaction 
2 and both transverse vector boson momenta. This is similar (but not identical) to single 
Drell-Yan production, where a correlation between the transverse polarization of a parton 
and its transverse momentum induces an azimuthal correlation between the momenta of 
the vector boson and its decay lepton [64]. 

We finally turn to the case where both vector bosons are produced from transversely 
polarized quarks. The corresponding contribution to the cross section is 




e 



%Z\q\—lZ 2 q% 



I 



d 2 y 



X I [C0S2((^1 - if 2 ) (K Sqi Sq 1 K Sq2 8q 2 + K'sq^K'sq^) 

- Sin2(</?1 - lf 2 ) (K'sqrffcKsqzSqi ~ Ksq^K'sq^)] 



+ [shl2((£i - if 2 ) (Ksq 1 Sq 1 Ks q2 5q 2 + K'sq^K'sq^) 
+ COs2((£i - (f 2 ) (K's qi Sq 1 K S q 2 Sq 2 ~ Ksq^K'sq^)] 



X y 2 M 2 (f&qi^dSq^ + SSq^fiqiS^) 



+ [C0S2(</?1 + (f 2 - 2(fy) (K S q 1 8q 1 K S q 2 Sq 2 ~ K5qi5qi K Sq2Sq 2 ) 
- Sin2(</?1 + (f 2 - 2tp y ) (K's qi Sq 1 Ks q2 Sq 2 + Ksq 1 8q 1 K's q2 Sq 2 )] 



y lvl \J SqiSq2 J SqiSq2 "Sq±8q2 "5<?i<5q2 / 



- [sin2(<£i + if 2 - 2ip y ) (Ks qi 8q 1 Ksq 2 Sq 2 - K8q 1 6q 1 K5q 2 8q 2 ) 
+ COS 2((pi + if2 — 2ify) (K'sq^KsqzSqz + Kgq^K'sq^)] 




+ {flavor interference} + {qq permutations} 



(3.4) 



-14- 



and depends on the azimuthal angles ipi, ip% and ip y in addition to the polar angles 0\ and 
02- The flavor interference and qq permutation terms are again obtained as in (3.2). 

The terms depending on ifi — <^ 2 describe a transverse correlation between the leptonic 
decay planes of the vector bosons. By contrast, the terms with ip\ + (f2 — 2ip y describe 
an azimuthal correlation between the lepton momenta and the direction between the hard 
interactions, which after integration over y, Z\ and 2 2 turns into an azimuthal correlation 
between the lepton momenta and the momenta of the two bosons. 

3.1 Cross section integrated over transverse boson momenta 

Integration over the transverse momenta of the two vector bosons yields cross sections ex- 
pressed in terms of collinear double parton distributions F(x±, X2, y)- Their spin structure 
is as in (2.9) to (2.11) but without pseudoscalar functions g, because one cannot construct 
a pseudoscalar with only one vector y. 

Upon integration over q l and q 2 , the cross section (3.1) for W pair production becomes 



d ° ww y k -k - 

X |(1 + COS^l) 2 (1 + COS# 2 ) 2 J d 2 y (f qiq2 + /AgiAgaXiW + /Ag 3 Ag 4 ) 
+ (1 + COS#i) 2 (1 - COS# 2 ) 2 J d 2 y (f qiq - 4 - fAq 1 Aq 4 )(fq 3 q 2 ~ fAq 3 Aq 2 ) 
+ (1 - COs6»i) 2 (1 + C0S6> 2 ) 2 J d 2 y (f qsq2 - fAq 3 Aq 2 )(f qi q 4 ~ fA qi Aq 4 ) 
+ (1 - COS#i) 2 (1 - C0S6» 2 ) 2 J d 2 y (fq m + f Aq 3 Aq 4 ){f qi q 2 + /AgiAg 2 ) 

+ {flavor interference} , (3-5) 

where the arguments of the distributions are f(x\,X2,y) and f(x\,X2,y)- In the general 
case we have a contribution 

tfcrW _ 1 

n 2 =i dxi dxi dQi c qiq2qMi 



[(1 + COS 2 e^Kg^ + 2 COS 0! K' qiq3 ] [(1 + COS 2 9 2 )K q2 g 4 + 2 cos 2 K' q2qi ] 
X [ d 2 y (f qi q 2 fq 3 q 4 + /a 9i Aq 2 fAq 3 Aq 4 ) 
+ [(1 + COS 2 0l)K qi Aq 3 + 2 COS 0i K' qiA q 3 ] [(1 + COS 2 2 )K q2 A q4 + 2 COS 2 K' q2 Aq 4 ] 
X J d 2 y (f qi q 2 fAq 3 Aqi + /a 9 i Ag 2 /j394 ) | 

+ {flavor interference} + {qq permutations} (3-6) 
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from unpolarized and longitudinally polarized partons. The contribution with transverse 
quark polarization in one of the two hard interactions now vanishes, 

**! =0. (3.7) 

lli=i dxi dxi dTti 

This is because integration of (3.3) over q l and q 2 sets z\ = z 2 = 0, after which the y 
integral gives zero due to the azimuthal dependence on ip y . By contrast, the contribution 
with transverse quark polarization in both hard interactions remains nonzero, 

— 2 sin 2 0i sin 2 2 



n.j=i dxi dxi dVLi C 



HQ2 



X | [COS 2(lfi - (f 2 ) (Ks qi 5q l Ks q2 8q 2 + K'sq^K'sg^) 

- sin2(v?i - <p 2 ) (Ks qi s qi K Sq2 s q2 ~ K SqiSqiK's q2 sq 2 )\ j d 2 y fs qi s q2 fs qi s q2 

+ {flavor interference} + {qq permutations} . (3-8) 

According to (2.11) the distribution fs qi 5 q2 describes the correlation between the directions 
of the transverse polarizations of two quarks in the proton. This correlation and its coun- 
terpart for antiquarks induce a correlation between the leptonic decay planes of the vector 
bosons, even if their transverse momenta are integrated over. Only if one integrates over 
the azimuthal angle of at least one of the leptons does the contribution from transverse 
quark polarization completely disappear from the cross section. 

The cross section of the double Drell-Yan process with two photons was calculated 
in [46], integrated over the transverse boson momenta and over the angles of the decay 
leptons. The expression in equation (9) of [46] agrees with our result (3.6) (up to the 
combinatorial factor 1/C, which was omitted in [46]). 

4 Transverse position dependence of distributions 

So far we have focused on spin correlations in DPDs and their consequences for the double 
Drell-Yan cross section. Even for unpolarized partons, however, there can be correlations 
between two partons in the proton, namely correlations affecting the dependence of DPDs 
on the transverse variables y, z\ and z 2 and the interplay between these variables and the 
longitudinal momentum fractions. 

In the present section we take a brief look at this issue by using a simple model in 
which the proton is described by a three-quark wave function. This is clearly too simple 
to describe the physics of small momentum fractions most relevant at the LHC, although 
it may actually be used for modeling quark DPDs at momentum fractions in the valence 
region. We proceed with this model in the spirit of an exploratory study. 

Our model ansatz for the three-quark light-cone wave function of the proton is 



*(»i,6i - b) = exp 



1 3 



4a 2 



(4.1) 
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where a is parameter of dimension length, b = x\bi+x 2 b 2 +x^b2, is the transverse position of 
the proton, and x\ + x 2 + X3 = 1. The corresponding wave function depending on transverse 
momenta is a Gaussian with exponent —a 2 k 2 /x i , which was long ago proposed in [65] 
and is often used for the phenomenology of valence dominated quantities, see e.g. [66] . The 
relation between the light-cone wave functions in transverse momentum and transverse 
position representation can be found in [67]. We do not specify the longitudinal part $(xi) 
of the wave function nor its spin-flavor dependence here, since the focus of our study is on 
the transverse variables. 

From the light-cone wave function (4.1) one can compute the contribution of the three- 
quark Fock state to the DPD of two quarks in the proton, in full analogy to the well-known 
case of single-parton distributions (discussed e.g. in [67]). Up to a factor depending only 
on the longitudinal momentum fractions Xi, the DPD is given by 



F(xi,Zi,y) oc exp 



1 f x-2 o™ „ - - , „ ,1 _ ^2 , 4xiX 2 n 



, xi(l - x\)z 1 - 2x x x 2 z x z 2 + x 2 {l - x 2 )z 2 H ■ y 

aa z I x\ + x 2 



d 2 b exp 



1 Xl+X 2 ( x\ 

— o -| y 

2a z 1 — x\ — x 2 V x\ + x 2 



(4.2) 



where b is the transverse position of the proton, averaged over the scattering amplitude 
and its conjugate as specified in [40]. The second line in (4.2) just gives an Xi dependent 
factor after integration over b. 

Inserting (4.2) into the cross section formula (2.1) and performing the integrals over 
all transverse positions, one obtains a cross section for double hard scattering that depends 
on the transverse boson momenta as 

exp< —a 2 



q 1 C n (x i ,x i ) + 2q 1 q 2 C 12 {x i ,x i ) + q 2 C 22 {xi,Xi) ^ j (4.3) 

with dimensionless functions Cj,- of the momentum fractions x±,x 2 and x±,x 2 , which are 
somewhat lengthy and will not be given here. The expression in square brackets is positive 
definite, so that the transverse momentum dependence has a Gaussian falloff at large 
transverse momenta. The coefficient C\ 2 describing the correlation between q l and q 2 is 
positive as well, so that one finds a preference for the two vector bosons to have opposite 
transverse momenta. We see that even with the simple wave function ansatz (4.1) the 
dependence of the cross section on the transverse momenta of the gauge bosons is not 
independent of their longitudinal momenta. 

An ansatz often made in phenomenology is to neglect correlations between partons 
and to write a collinear DPD as the convolution of two single-parton distributions that 
depend on the momentum fraction and the transverse position of the parton. This ansatz 
can be extended to include Zi dependent DPDs and then reads 

F(xi,Zi,y)^ j d 2 b f{x 2 ,z 2 \b+\xizi) f{xi,zi]b + y -\x 2 z 2 ) . (4.4) 

The second argument of the single-parton distribution / is Fourier conjugate to the trans- 
verse quark momentum and the third argument gives the transverse position of the proton 
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with respect to the quark, both averaged over the scattering amplitude and its conjugate. 
The shift of this argument by ^x\Z\ or — \x 2 z 2 is a consequence of Lorentz invariance 
as explained in [39]. Evaluating the single-parton distributions for the light-cone wave 
function (4.1) one obtains 



F(xi,Zi,y) oc J d 2 b exp 



oc cxp 



x exp 
1 



8^T^7{ (1 -^ + (26 + 2! ™> 2 } 



xi(l - xx)z\ + x 2 (l - x 2 )z\ 



+ 



X\X 2 



- X 2 ) + (1 - X\)x 2 



2y - x\Z\ - x 2 z 2 



(4.5) 



for the transverse dependence of the DPD. This is visibly different from the result (4.2) 
of the direct calculation. Although the ansatz (4.4) involves the convolution of two single- 
parton distributions, thus suggesting that the two partons are distributed independently, 
it induces correlations between transverse and longitudinal variables in F(xi,Zi,y). 

Inserting the form (4.5) into the cross section formula one obtains again a Gaussian 
behavior as in (4.3), but with different coefficients Cij. In particular, the sign of C\ 2 is then 
equal to the sign of {x\ — X\){x 2 — x 2 ), so that depending on the longitudinal momentum 
fractions the transverse boson momenta tend to be in the same hemisphere or in opposite 
ones. This difference in qualitative behavior shows that the ansatz (4.4) must be used with 
great care when one is interested in correlation effects. 

Setting z\ = z 2 = in (4.2) and (4.5) gives collinear DPDs with a Gaussian depen- 
dence on y. The Gaussian width depends on x± and x 2 and differs in the two cases, 



F(xi,y) 
F(xi,y] 



(4.2) 



oc exp 



1 



X\X 2 n 

y 



(4.5) 



oc exp 



2a 2 x\ + x 2 

1 XiX 2 



2a 2 x\ + x 2 — 2x\x 2 



y~ 



(4.6) 



We see that, within our model, the ansatz (4.4) does not reproduce the interplay between 
y and the momentum fractions. It does, however, provide a valid approximation unless x\ 
and x 2 are both rather large. 



5 Summary 

Multiple hard interactions in pp collisions can yield substantial contributions to the pro- 
duction of final states with high multiplicity in parts of phase space. In this paper we have 
shown how spin correlations between two partons in the proton affect the rate and the 
angular distribution of the final state in the production of four leptons via two electroweak 
gauge bosons. We considered both the case where the transverse momenta of the bosons 
are small (using transverse-momentum dependent factorization) and the case where they 
are integrated over (using collinear factorization). 
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We find that longitudinal spin correlations between the quarks or antiquarks in the pro- 
ton affect the rate of double parton scattering, and in the presence of axial- vector currents 
also the polar distribution of the produced leptons. Correlations involving transversely po- 
larized quarks or antiquarks induce azimuthal correlations between the final state leptons. 
A part of these correlations persists if the transverse momenta of the gauge bosons are 
integrated over. Having two "independent" hard interactions in double parton scattering 
does hence not imply that the final states produced by the two interactions are independent 
of each other. 

How large parton spin correlations inside a proton actually are remains an open ques- 
tion that deserves further study. This also holds for possible correlations between the 
transverse distribution and the longitudinal momentum fractions of the partons. We find 
several such correlations in a simple model with a three-quark wave function. Within this 
model we also find that the often used ansatz to represent double parton distributions as 
convolutions of single-parton distributions is inadequate to describe details of the kinematic 
dependence in double parton scattering. 

Double Drell-Yan production involves a particularly simple hard-scattering subprocess 
but nevertheless exhibits a rich pattern of angular effects induced by parton spin correla- 
tions. It is natural to expect that other processes, in particular those involving multijets, 
will share this feature. We note that the cross section dependence on angles between the 
final-state particles implies a dependence on the invariant mass of particle pairs, which is 
an important quantity in searches for new physics. An estimate of the possible size of such 
effects would therefore be of great value. 
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A Coupling Factors 

In this appendix we list the coupling factors K and K 1 appearing in the double Drell-Yan 
cross section. Further relations between these factors are given in section 2.3. 

A.l Charged vector bosons 

For W + production one has 



a 



2 IT/ |2 



W Q 



K ™> ~ AN C (2 sin 9 w r (Qf - m^f + m^ w ' (e « e * " 1} ' (AJ) 
and for W~ production 



q 2 \v mi \ 2 Qf 

4iV c (2 sin^,) 4 (Qi ~ m 2 w ) 2 + 



Kqiqj a at m „:L a \4 I r\2 __2 \2 , „_2 r2 ' ( e <?i e <?i — 1 (A-2) 



Here N c = 3 is the number of colors, V qiqj a CKM matrix element, 8 W the weak mixing 
angle, a the electromagnetic fine structure constant, and e qi the charge of quark qi in units 
of the positron charge. 
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A. 2 Neutral vector bosons 

For a lepton pair £ + £~ produced via a 7*, Z 01 their interference, one has coupling factors 

2 f 6 6 

K im = ^rj^)^ - MQi) 9i (fiqi9q + %4) - iB(Qi)gY(e qi gl - e q .g^.) 

+ Cm [{gjf + (4) 2 ] (44 + 44)} , 
K' m = ^-\-A(Qi)9lKgf. + e 9 ,4) - iB(Q<) gf(e qi gf. - e qj g£) 

+ C{Q i )2gJgt{glg^ j +g v qj g^ 



# ft A$ = ~A(Qi) gj{e. qi g qi + e 9j 4) _ iB {Qi) 4( e <7i4 - e Qj g 

+ C{Q t ) [(gYf + (4) 2 ] (44 + 44) 



K'^Aqj = jjjrl ~A(Qi) gi(e qi g v qj + e qj g^.) - iB{Qi) gf{e. qi g v q . - e % g v qi ) 

+ C(Q0 2^4(44 + 44)1 ( A - 3 ) 



and 



2 r g g 

= -^l 1 - MQi) 4( e 9»4 + e ^4) - t-WO 4(^4 ~ %4) 

+ c(g 4 ) [(4) 2 + (4) 2 ] (44 - 44)} > 
2 ( 

K' 8qiSqj = ^\-B{Qi)gY{e q .g^ + e^gfy +iA(Qi) gj (e^gf. - e^g*) 

- iCm [(gjf + {gff] (44 - 44)1 . (A.4) 



Here we have used the conventional vector and axial fermion couplings to the Z boson, 

gj = I 3 f -2e f sm 2 e w , gf = ij , (A.5) 

where J| is the third component of the weak isospin of the left handed fermion / and ej 
its charge in units of positron charge. Since we do not consider Z decays to neutrinos, £ is 
always a negatively charged lepton. We have furthermore used the abbreviations 

A{Qi) — — 2 0/1 77-i2 2 \9 I 2~FT ' B(Qi) 



sin 2 26 w (Q 2 - m|) 2 + m| T% ' ^ sin 2 20^ (Q 2 - m|) 2 + m|r| ' 
1 Q 2 

° m = sln^ (Qf-m|) 2 + m|r| ' (A ' 6) 

For the usual hard-scattering cross sections one has equal flavors qi = qj in the above 
coupling factors and finds that their imaginary parts are zero. This is not the case for the 
coupling factors describing flavor interference, where qi 7^ qj. 
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